o 

(N 

ON 



Nonlinear thermal transport and negative differential thermal conductance in 
graphene nanoribbons 

Jiuning Hu, 12 'E3 Yan Wang, 3 Ajit Vallabhaneni, 3 Xiulin Ruan, 3 ' 2 and Yong P. Chen 4 ' 2 ' 1 ^ 

^School of Electrical and Computer Engineering, Purdue University, West Lafayette, Indiana 47907, 
USA 

2 t Birck Nanotechnology Center, Purdue University, West Lafayette, Indiana 47907, 
USA 

^School of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47907, 
USA 

^Department of Physics, Purdue University, West Lafayette, Indiana 47907, 
USA 

We employ classical molecular dynamics to study the nonlinear thermal transport in graphene nanoribbons 
(GNRs). For GNRs under large temperature biases beyond linear response regime, we have observed the onset 
of negative differential thermal conductance (NDTC). NDTC is tunable by varying the manner of applying 
the temperature biases. NDTC is reduced and eventually disappears when the length of the GNR increases. 
We have also observed NDTC in triangular GNRs, where NDTC exists only when the heat current is from 
the narrower to the wider end. These effects may be useful in nanoscalc thermal managements and thermal 
signal processing utilizing GNRs. 



o 

(N 
> 

oo 
o> 
o 

o 



X 



GraphcnCfi^ an atomic monolayer of graphite, has 
emerged as one of the most interesting materials in 
condensed matter physics and nanotechnology. Besides 
its unusual electronic properties^ graphene also has 
unique thermal properties, e.g., high thermal conduc- 
tivities (~600-5000 W/m-K)£~— Graphene nanoribbons 
(GNRs) are promising in many applications, such as their 
electronic band-gap tunability^ and edge chirality depen- 
dent thermal transport^ So far, little attention has been 
paid to nonlinear thermal transport in GNRs, though 
these nonlinear effects have been explored in ideal atomic 
chainsjiSr— molecular junctions^ and quantum dotsJ£ 
Here, we demonstrate negative differential thermal con- 
ductance (NDTC) in GNRs. Analogous to the electronic 
counterpart fii NDTC is a useful ingredient in developing 
GNR-based thermal management and signal manipula- 
tion devices, such as the thermal amplifiers^ and ther- 
mal logic gates J£ 

We study the thermal transport in GNRs using clas- 
sical molecular dynamics (MD) simulations. The many- 
body empirical Brenner potential is employed to de- 
scribe the carbon-carbon interactions. This method have 
been applied in many graphene-based systems ^i^2r— The 
structures of GNRs in this study are shown in the in- 
set (rectangular GNR) of Fig. Q] and the inset (triangu- 
lar GNR) of Fig. [3] The atoms denoted by squares are 
fixed in position, while those denoted by left- and right- 
pointing triangles are placed in two Nose- Hoove r 23 ' 24 
thermostats at different temperatures Tl and T R , respec- 
tively. The equations of motion for atoms without posi- 
tion being fixed are: 
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where Pi is the momentum of the i-th atom, Fj is the 
total force acting on the z-th atom, and ji is the Nose- 
Hoover dynamic parameter. For the atoms denoted by 
circles, 7; = 0, and it recovers the NVE (constant number 
of atoms, volume, and energy) ensemble. For the atoms 
in the left and right thermostats, 7, obeys the equation 
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where r is the thermostat relaxation time, N 
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number of atoms in the thermostat, fcs is the Boltzmann 
constant and m is the mass of the carbon atom. More 
details on our numerical calculation method can be found 
elsewhere 

First, we study the thermal transport in a rectangu- 
lar GNR with armchair top and bottom edges shown in 
the inset of Fig. [1] (we have obtained qualitatively similar 
conclusions for GNRs with zigzag edges) . Since the GNR 
is symmetrical, we only consider Tl < T R and define the 
temperature difference AT = T R — Tl ■ The temperature 
T R is kept as a constant. As we can see from both curves 
in Fig. [U for small temperature difference (e.g., AT < 60 
K for T R = 300 K and AT < 150 K for T R = 600 K), 
the thermal current increases approximately linearly as 
AT increases, as expected from Fourier's law. Interest- 
ingly, for some range of higher AT, the thermal current 
decreases as AT increases (the dashed boxes in Fig. [J), 
indicating the onset of NDTC. It is a reasonable approxi- 
mation to consider thermal current as proportional to the 
product of thermal conductivity k of the GNR and AT. 
Our previous study^ has shown that k increases with the 
average temperature T = (T L + T R )/2 = T R - AT/2. 
We have plotted T (labeled at the right vertical axis and 
indicated by right-pointing arrows for Fig. 1-3 and in the 
subplot of Fig. SJb)) as a function of AT in all figures. 
Since T decreases with AT, k decreases with increasing 
AT. The resulting trend of the thermal current as a func- 
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FIG. 1. Thermal current (left vertical axis) and average tem- 
perature (right vertical axis) vs. temperature difference AT. 
The dashed boxes highlight NDTC. The inset shows the struc- 
ture of the GNR (~ 1.5 nm x 6 nm). ■ denotes fixed bound- 
ary atoms. -4 (►) denotes atoms in the left (right) thermostat. 
• denotes the remain atoms in the bulk. 



FIG. 3. Thermal current (left vertical axis) and average tem- 
perature (right vertical axis) vs. temperature difference AT 
in triangular GNR shown in the inset. The labels for the 
GNR structure have the same meaning as that in the inset in 
Fig. [I] The dashed box highlights NDTC. 
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FIG. 2. Thermal current (left vertical axis) and average tem- 
perature (right vertical axis) vs. temperature difference AT 
in GNRs with the similar structure as the GNR in the in- 
set of Fig. [T] except for different lengths. In all these plots, 
T R = 300 K and T L is varied from 300 K to 30 K. 



tion of AT is thus a competition between decreasing k 
and increasing AT. In the AT range displaying NDTC, 
the decrease of k with AT dominates. We have found 
that there is no NDTC (shown in Fig. 0]) if T L is larger 
than the constant Tr, i.e., if T increases with AT (thus 
without the above competition). Note that for large AT 
beyond linear response, strictly speaking thermal con- 
ductivity is not well defined. Thus, in the above expla- 
nation, k is considered to be an effective, average ther- 
mal conductivity. Similar arguments have been applied 
in analysing thermal transport in ID atomic chains.— 

Second, we study the length dependence of NDTC in 
GNRs. For all three GNRs of different lengths in Fig. d 



T R = 300 K while T L is varied from T R to 30 K. As the 
GNR length is increased, the AT value for the onset of 
NDTC increases and the AT range where NDTC exists 
shrinks. We thus suggest that NDTC will eventually dis- 
appear if the length of GNR exceeds some critical value. 
We have verified this using LAMMPS package^ and ve- 
locity scaling^ MD, and found no NDTC in a 50 nm long 
GNR with similar structure as that studied in Fig. Q] 

Besides these nonlinear effects in symmetrical GNRs, 
we also explore the possibility of NDTC in an asymmet- 
rical triangular GNR, shown in the inset of Fig. [3] Our 
previous study has pointed out that thermal rectification 
exists in this asymmetrical GNRi^ As we see from Fig. EH 
here the nonlinear thermal transport is also direction- 
dependent. NDTC appears when the temperature of the 
narrower end is held at Tl = 300 K and the tempera- 
ture Tr of the wider end is varied from 300 K to 30 K 
(solid line in Fig. [3]). However, there is no NDTC when 
the values of Tl and Tr are interchanged (dashed line 
in Fig. [3]). This provides another possibility to control 
the nonlinear thermal transport and NDTC in GNRs by 
engineering the shape of GNRs. 

In general, the way to tunc the thermal current in the 
two-terminal thermal devices is very different from that 
in any two-terminal electronic devices. In the latter case, 
only the voltage difference matters. However, in thermal 
devices, the average temperature T is as important as 
the temperature difference AT in controlling the ther- 
mal current. For example, consider T = aAT + To with 
constants a and To, and we have Tl — (a — ^)AT + To 
and Tr — (a + §)AT + T . The thermal currents and 
average temperature T as a function of AT are plotted 
in Fig. |4] for the rectangular GNR shown in the inset in 
Fig. [[J where T = 300 K and a is tuned from -0.5 to 0.5 
(indicated by the dashed curved arrow in Fig. 3]). The 
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FIG. 4. Thermal current (a) and average temperature (b) vs. 
temperature difference AT for different values of a for the 
GNR shown in the inset of Fig. Q] Note that a = 0.5 (-0.5) 
corresponds to T L ( R ) fixed at 300 K while Tml) is varied. 



solid curve in Fig. [T] corresponds to a = —0.5. For small 
temperature difference in the linear response regime, the 
slope of thermal current vs. AT is independent of a. 
In the nonlinear response regime (large AT) , the system 
transitions from a regime with NDTC to a regime with- 
out NDTC when a is tuned from negative to positive 
values. We can see a strong correlation between the the 
trend of the thermal current and that of the average tem- 
perature for different values of a in the range of AT from 
100 K to 250 K where NDTC occurs for negative a. For 
negative a, since T decreases with AT, the effective k 
decreases with AT, and the occurrence of NDTC can be 
similarly explained as that for Fig. Q] 

There are two independent parameters to control 
the thermal transport in two-terminal devices, cither 
(Tl,Tr) or (AT, T). Two-terminal thermal devices are 
actually analogous to three-terminal electronic devices. 
In the language of electronic transport of field effect tran- 
sistors (FETs) , AT plays the role of the drain-source volt- 
age difference in FETs, while a plays the role of the gate 
voltage. Fig. @] shows the ability to realize the FET-like 
behaviour in GNRs. 

In summary, we have studied the nonlinear thermal 
transport in rectangular and triangular GNRs under 
large temperature biases. We find that in short (~ 6 nm) 
rectangular GNRs the NDTC exists in a certain range of 



applied temperature difference. As the length of the rect- 
angular GNR increases, NDTC gradually weakens. In 
triangular GNRs, NDTC only exists in the thermal cur- 
rent direction from the narrower to the wider end. The 
ability to tunc and control NDTC by temperature param- 
eters and GNR shapes provides potential ways to manage 
heat and manipulate thermal signals at the nanoscale. 
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